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Abstract 

We point out an elementary thermodynamics fact that whenever the specific heat of a 
system is negative, the speed of sound in such a media is imaginary. The latter obser- 
vation presents a proof of Gubser-Mitra conjecture on the relation between dynamical 
and thermodynamic instabilities for gravitational backgrounds with a translationary 
invariant horizon, provided such geometries can be interpreted as holographic duals 
to finite temperature gauge theories. It further identifies a tachyonic mode of the 
Gubser-Mitra instability (the lowest quasinormal mode of the corresponding horizon 
geometry) as a holographic dual to a sound wave in a dual gauge theory. As a specific 
example, we study sound wave propagation in Little String Theory (LST) compacti- 
fied on a two-sphere. We find that at high energies (for temperatures close to the LST 
Hagedorn temperature) the speed of sound is purely imaginary. This implies that the 
lowest quasinormal mode of the finite temperature Maldacena-Nunez background is 
tachyonic. 



July 2005 



1 Introduction 



Gauge theory/string theory correspondence [1,2] allows to rephrase complicated ques- 
tions arising in string theory (gravitational) backgrounds into more intuitive field the- 
oretical language. In certain cases this dual formulation of the problem allows for a 
simple solution. In this paper we discuss a problem of just this type — the Gubser- 
Mitra conjecture [3]. 

The Gubser-Mitra conjecture^ states that gravitational backgrounds with a transla- 
tionally invariant horizon (the 'black brane' solutions) develop an instability (a tachy- 
onic mode) precisely whenever the specific heat of the black brane geometry becomes 
negative. In the framework of Maldacena duality gravitational backgrounds with trans- 
lationary invariant horizon arise as a dual description of strongly coupled generalized^ 
gauge theories at finite temperature. In what follows we consider only black brane ge- 
ometries of the latter type. The holographic dual of the Gubser-Mitra conjecture then 
implies that a finite temperature gauge theory with a negative specific heat must have 
a dynamical instability. It is straightforward to explicitly identify such an instability^. 
Indeed, consider a thermodynamic system (without chemical potential) at fixed volume 
and temperature T. Such a media propagates small momentum (g <^ T) sound waves 
with a dispersion relation 

u:{q)=Vsq + o(^^^ , (1) 
where the speed of sound Vs is determined from the equation of state as 

"-'as' 

with P, S being the pressure and the energy density correspondingly. At zero chemical 
potential, the free energy density JF is JF = —P. Noting that at a fixed volume V 

related to the entropy density S, and using the specific heat cy definition 

^Interesting observations on Gubser-Mitra conjecture were presented in [4,5]. 

^In some examples [6] the Maldacena duals of string theory backgrounds have a four dimensional 

gauge theory interpretation only at low energies. 

After completing this paper we learned that equivalent arguments were presented in [7]. 
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we find from Q 

^ {dP/dT)y ^ S_ 

' {d£/dT)v cv' ^ ' 

Since^ 5* > 0, eq. (0) implies tliat in a media witli cy < tlie speed of sound is purely 
imaginary. Thus, the amplitude of a fixed momentum sound wave would increase 
exponentially with time, reflecting the dynamical instability of the system. Finally, 
as the sound wave in hot gauge theory plasma is holographically dual to the lowest 
quasinormal mode of the dual black brane geometry [8], we immediately conclude that 
it is this mode that realizes tachyonic instability of the black brane backgrounds with 
a negative specific heat. 

In the rest of this paper we explicitly demonstrate above observation in the con- 
text of Little String Theory [9, 10] compactified on S"^. The extremal supergravity 
solution representing a large number of type IIB string theory NS5 branes wrapped 
on a two-cycle of the resolved conifold was discussed in [11,12,6] (MN). In the in- 
frared MN geometry provides a holographic dual to A/" = 1 SU{N) supersymmetric 
Yang-Mills theory in the planar limit. This supergravity solution encodes interesting 
non-perturbative phenomena of the four- dimensional asymptotically free gauge theo- 
ries: confinement, chiral symmetry breaking, etc. The nonextremal deformation of the 
MN background (holographically dual to 5*^ compactified LST at finite temperature) 
was previously discussed in [13-15]. We compute the dispersion relation for the lowest 
quasinormal mode in the finite temperature MN geometry, as described by (P). We 
find that at high energy (at temperatures close to the LST Hagedorn temperature) the 
speed of sound is purely imaginary. We then reproduce the speed of sound directly from 
the equation of state^ computation (0). This provides a highly nontrivial consistency 
check on our analysis. 

One of the motivations to consider finite temperature deformations of the Maldacena- 
Nunez background [6] was to study finite temperature phase transitions (confine- 
ment /deconfinement, chiral symmetry breaking/restoration) in a QCD-like theory. Un- 
fortunately, we demonstrate here that such deformations are tachyonic^ for the MN so- 
*in the gravitational dual it is proportional to the area of the horizon 

'^Though both Ref. [14, 15] propose the (different) equation of state for the finite temperature MN 
geometry, neither of them reproduces the correct speed of sound. We comment on this issue later in 
the paper. 

^Strictly speaking, we establish this only at high energy. We believe that this is true at all energy 
scales where supergravity approximation to the full string theory is reliable. 
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lution. The only currently known string theory background where thermal phase tran- 
sitions in a Q CD-like theory can be studied quantitatively is the Klebanov-Strassler 
geometry [16-21]. 

The paper is organized as follows. In the next section we derive effective five 
dimensional gauged supergravity action describing finite temperature deformations of 
the MN background at temperatures above chiral symmetry breaking. In section 3 
we construct analytic solution of the nonextremal geometry holographically dual to 
S"^ compactified LST, in the limit of large 5*^ radius at the horizon. Wc argue that 
this parameter regime corresponds to the high energy phase of the finite temperature 
MN solution. We motivate a new equation of state for the system and explain the 
difference with the previous proposals [14,15]. In section 4, following [8,22], we compute 
the dispersion relation for the sound wave in finite temperature MN background. We 
conclude in section 5. 



2 Effective action for LST on 

Extremal and nonextremal solutions of the type IIB supergravity backgrounds holo- 
graphically dual to S"^ compactified LST with A/" = 1 supersymmctry preserving twist 
(at the extrcmality) were constructed in [11,12,6,13,15]. In this section we present 
an effective five-dimensional supergravity action describing such solutions (and pertur- 
bations about them) assuming the unbroken U{1) symmetry of the MN background 
(corresponding to the U{1)r symmetry of the dual low-energy A/" = 1 supersymmetric 
Yang-Mills theory). This effective action provides a consistent five-dimensional trun- 
cation of the full type IIB supergravity geometry of a large number of NS5 branes 
wrapped on a two-cycle of the resolved conifold^. 

Consider the following type IIB supergravity ansatz. The string frame metric is 

^^strtng _ ^ ^^^rf-tJ- (^rf-f _|_ (^^g^y^ _|_ _|_ ^^2 _|_ cOS 6*1 , (6) 

where uj — a;(x^), {^^i, 0i} are coordinates of the round 5*2, {6'2, ^2} are coordinates of 
the round 5*2, n is a constant. Additionally we have a 3- form flux 

H2, = i^i^ + ^ cos 9i d^H^ A (sin Oi dOi A d^i - sin 6*2 d02 A #2) , (7) 

^Such cfFcctivc action is a direct analogy of the 'chirally symmetric' action for the Klebanov-Tseythn 
geometry constructed in [21]. 
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and a dilaton $ = <l>(x^). Evaluating type IIB (Einstein frame) supergravity action 

levrGio Jmio V 2 2 2 

- ^Fs A ^Fs - AHsAFs 

on the ansatz we obtain five dimensional effective action dual to LST compact- 

ified on S'^ 

IbvrGs y^vig I 
+ + 2cu"^ - n^cj^^l , 

(9) 

where the greek indexes run 1...5. is the 5d Ricci scalar computed with the metric 

dsl = g^y{x)dx^dx'' , (10) 
G5 is a 5d effective gravitational constant 

We find it convenient to rewrite the action © in five- dimensional Einstein frame. 
The latter is achieved with the following rescaling 

g,u^9,u^e-""'u'/' g,,. (12) 

Further introducing 

a = Ina; , (13) 
the five dimensional effective action becomes 



If ( 10 o 4 9 8 

/ ^^^^^ \ ^^-^ (^") - ^ (^'^) + 
l67rG5 Jj^^ 1^ 3 3 3 

„ 4. 4 4. 10 „ 4. 16 

Prom the effective action ()14|) we obtain the following equations of motion 

4^ 10 ^ 4, 16 



(14) 



1 4^4 1 4^ 16 

=□$ + -n-2e3*-3- + -n^es*--" , (15) 
10 4 4 1 



where we denote 



^ = ^-2e3*-3- + 2e3*-T 




(16) 



3 Finite temperature MN geometry: analytic solutions and 
the equation of state 

We begin with deriving equations of motion describing finite temperature deformation 
of the chirally symmetric MN solution [13, 15]. The regular horizon solutions of these 
equations are characterized by two independent parameters: the value of a dilaton and 
the 5*2 (see (jH))) radius at the horizon. These parameters are related to the temperature 
and the strong coupling scale^ of the dual SYM theory (alternatively the Kaluza-Klein 
scale of the compactified LST). We analytically construct solutions for the nonextremal 
deformations perturbatively in inverse 5*2 radius at the horizon, which corresponds to 
a high energy (near Hagedorn ) regime of the compactified LST [14]. We argue that 
keeping the strong coupling scale of the dual gauge theory fixed corresponds to a very 
specific dependence of the horizon value of a dilaton on an 5*2 horizon radius^. We 
finally derive high-energy equation of state for the compactified LST and discuss the 
discrepancy with the previous proposals [14,15]. 

3.1 Background equations of motion for the finite temperature MN solu- 
tion 

Nonextremal deformations of the MN geometry are described by the following ansatz 
of the five-dimensional effective action (jH)) 



*^The fact that only two parameters specify regular horizon solution implies that finite temperature 
gauge theory dual to nonextremal MN background indeed does not have a chemical potential, and 
thus falls in a class of systems discussed in section 1. 

^This differs from the assumption in [14] where the high temperature thermodynamics of the MN 
solution was deduced assuming constant horizon value of the dilaton. 




(17) 
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where / = f{r), also we have $ = $(r), a = ^\ng + Inn = ^lng{r) + Inn. In this 
case equations of motion (fT3j) become 

0=r + [\ng-2^]'f', 

0=g"+[\nf-2^]'g' + ^, ^^g^ 
0=<l>" + [ln/(7-2$r$' + i-^, 

8g 

where prime denotes derivative with respect to r. Additionally we have a first order 
constraint 

-.2 



0=[\ng]' 



1 1 
In 5(8/2 - $ 



$'[$-ln/r + ^^^. (19) 



We specifically choose such a parametrization that all n-dependence drops out in (fTR|l . 
()19|). In this parametrization the extremal chirally symmetric solution of MN [6] takes 
a very simple form 

f = l, g = r, $ = 00 - + ^Inr. (20) 

Notice from ()2()j) that MN asymptotics (as r — > +00) can be summarized in a coordinate 
independent way as 

<t> = -^g + 0{lng). (21) 

We propose that enforcing the asymptotic (PT|) even off the extremality is equivalent 
to keeping the Kaluza-Klein scale of the compactified LST fixed. We will show later 
that this implies a particular dependence of the horizon value of the dilaton on the 5*2 
radius, also evaluated at the horizon. 

Nonextremal solutions of ()18|1 are best analyzed using a new radial gauge 

x = f. (22) 

This gauge choice is nonsingular because /(r) is a smooth, monotonic function^" of r 
in (fT7j) . Furthermore, it reduces the number of independent functions describing the 
nonextremal MN background to two: g{x), $(x). With the radial coordinate as in ()22|) . 
the black brane horizon is at x = and the boundary is at x = 1. The background 
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We explicitly verified this fact. 
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equations of motion (fTHj) take the form 

^=9 / o , ^ TT [9 > + , o , o 7^ ^ 



^(^2 + 2^-1) a;(^2 + 2(7-1) a;(^2 + 2(7-1) ' 

0=$" + ^i^!±lL .^/^2 2x((7^ + 1)(7^ + g'- 2g' + 3g ^, jg' + l)g'ixg' + Ag) 
g2 + 2g-l^ ' xg{g^ + 2g-l) Axg\g^ + 2g-l) ' 

(23) 

where all derivatives are now with respect to x. Regularity at the horizon implies the 
following perturbative expansion as x 0+ [13] 

oo oo 

^7 = «o + I] «n , $ = (/)o + ^ p„ x^" , (24) 

n=l n=l 

where na^ and 0o are the radius of the 5*2 and the value of the dilaton at the horizon 
respectively. The first couple terms of the perturbative expansion take values 

_ 4ag(l - al) Aal{al + ?,al-l){al-l) , 



al + 2al + l 2 ^^^^ 

Notice that pi is not fixed; in fact, it is easy to see that a„ oc p", p„ oc p". The latter is 
related to the exact scaling symmetry of ()23j) : x \x with g, <I> kept invariant^^. While 
ao and (j)Q are independent physical parameters of the nonextremal solution (related 
to the temperature and the strong coupling scale of the gauge theory), pi is not. As 
explained in [13], for a fixed values of {ao, 0o}) the parameter pi must be fixed in such 
a way that the boundary of the nonextremal geometry is indeed at x = 1. This is 
equivalent to requiring that g +oo as x ^ 1_. 

As shown in [13,15] the global solution of (j^^ (with a regular horizon) is crucially 
sensitive to the value of oq. For oq > 1 the solution is singularity-free in the bulk, and 
approach asymptotically the MN solution (j2ip . For < 1 the solution has a naked 
singularity in the bulk. Finally, for ao = 1, the exact solution takes form 

(7 = 1, $ = 0o + ^ln(l-x2). (26) 

Clearly, does not asymptote to MN solution near the boundary. Quite interest- 
ingly, the ao = 1 nonextremal background has an exact string world-sheet descrip- 
tion [15]. 

^"'^This scaling symmetry is broken by the boundary, located at a: = 1. 
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3.2 Analytic nonextremal solutions of eq. for ^ 1 

In previous section we derived equations of motion describing nonextremal deforma- 
tions of the chirally symmetric MN background ((2S1)- The globally regular solutions are 
characterized by the horizon boundary condition oq > 1. Geometrically, na^ is the ra- 
dius of the 5*2 (see ©) of the string frame ten-dimensional geometry. Whenever ao > 1, 
the string frame size of 5*2 monotonically increases as one approaches the boundary, 
diverging at the boundary. This 5*2 is a two-cycle of the resolved conifold, that is being 
wrapped by the NS5 branes. Intuitively, the larger the ao, the more locally fiat the NS5 
branes are. Thus we expect that the +oo limit of the nonextremal MN solution 

to describe the non-extremal type IIB string theory NS5 branes. The thermodynamics 
of this system is Hagedorn [23] . We are unable to solve analytically for arbitrary 
values of oq. On the other hand, it is possible to solve these equations perturbatively 
in Oq ^. Such perturbative solution describes corrections to the Hagedorn thermody- 
namics of NS5 branes wrapped on a large two-cycle. It is this regime where we test 
the proof of the Gubser-Mitra conjecture presented in section 1. 
For large values of ao we seek solution to in the form 

9=% + 2^ 



n=0 

oo 







(27) 

^ J, , 'Pnix) 



n=0 



It is straightforward to substitute ansatz (jTfj) into (j2Sl), and solve resulting ODE's 
iteratively in n. Given a horizon boundary condition 
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4, $ 



= 00 , (2^ 

horizon 



horizon 

and the MN asymptotic, {Qn{x)^Vn{x)} are uniquely determined. In what follows we 
will need an explicit solution up to order n = 2. We find 

Go = -21n(l-x2), 



2 (30) 
Gi = 21n(l -x^) +4 dilog(2;2) + 81na; ln(l - x^) - -7r^ 
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V2 =2 ln(l - x^) - 111^(1 - x^) - 2 dilog(l - x^) , 

^2 = -61n(l-a;^) -81nx ln^(l - a;^) - 8 ln(l - x^) polylog(2, 1 - x^) (31) 
+ 8 polylog(3, l-x^) + 12 dilog(l -x^) + 2 - x^) - 8C(3) • 

From (j27j), (pUjl - lj^ it is obvious that the hmit +oo does not commute with 

the hmit x —>■ Thus, a correct order of hmits must be specified. In what follows 
we always assume that the oq limit is taken first. The motivation for such a choice 
comes from the fact that as +oo the temperature of the nonextremal MN solution 
(see dnH)) approaches fiat NS5 brane Hagedorn temperature. Moreover, such a scaling 
implies that in the strict qq +oo limit only n = terms of ()27|) must be kept. In 
this case 

g =al + Qo-^ +00 , 

1 (32) 
$=0o + Po = </,o + -ln(l-x2), 

which is precisely the geometry of the nonextremal fiat NS5 branes [23]. Finally, such 
a scaling produces sound wave dispersion relation (see section 4) in agreement with 
that of the nonextremal fiat NS5 brane system [30]. Thus, it is natural to identify the 
proposed regime with the near-Hagedorn regime of NS5 branes wrapped on a large two- 
cycle. It is conceivable that some other order of limits would produce a system with a 
different thermodynamic (and hydrodynamic) description; we were not able to find a 
coherent description other than what we present here. It would be very interesting to 
find such a description (if it exists). 

We conclude with section by noting^^ that in the prescribed order of limits the MN 
asymptotic (PT|) is indeed satisfied provided we identify 

00 = 00 - ^flo + ^(In «o) , (33) 
where 0o is independent of oq. 

3.3 Large oq thermodynamics of the finite temperature MN solution 

In this section we study the thermodynamics of the LST compactified on a large two- 
cycle of the resolved conifold. 
i^This follows from 
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Given the metric ansatz ()17|). and the regular horizon boundary conditions 
we find the Hawking temperature of the black brane solution by identifying its inverse 
with the periodicity of the Euclidean time direction 

1/2 

4?™ \ 2pi 



^=7^f-7^(l + «o^)) • (34) 



Next we compute the Bekenstein-Hawking entropy density of the geometry (fT7j) . We 
find the three-dimensional area of the horizon ^3 of the black brane to be 

^3 = n^age-'^o^s , (35) 

where V3 is the three-dimensional volume. The entropy density of the black brane is 

5 = oc n^e-'^^ . (36) 
In the limit Oq ^ 1 we can use perturbative solution ^I7\i . ((221)- (EU to determine 

Pi = -l + a^' + 0{a,'). (37) 
From (jnil), (jTTjl we see that in the limit ao — >■ +00 [14] 

T = Th{1 + af + O(ao ^)) , (38) 

where Th = 1/(4?™) is the Hagedorn temperature of the six- dimensional LST. 

To determine equation of state we need to compute the energy density (and the 
pressure) of the gravitational background (fTTj) . The correct way to do this is to study 
holographic renormalization [24-28] of the MN solution, much like equivalent problem 
recently solved for the cascading gauge theories [21]. In all cases we are aware of, 
the properly implemented holographic renormalization in the context of gauge theory/ 
string theory correspondence produces the entropy density and the pressure automat- 
ically satisfying the standard thermodynamic relations^^ 

P = -T, dS = TdS. (39) 

Lacking the holographic renormalization for the MN background, one could still deduce 
the remaining thermodynamic quantities (besides the entropy density) by enforcing 
(I39p . The subtlety of the latter prescription in the context of MN geometry is that 
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For asymptotically locally AdS backgrounds this was rigorously proven in [29] . 
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a priori, one does not know whether 0o depends on oq once the physical properties 
of the gauge theory (in our case the strong couphng scale) is kept fixed, and if so, 
what is the precise relation. Identical problem arises in the study of thermodynamics 
of the cascading gauge theories [18-21]. There, the entropy density also depends on 
two parameters: the five-form fiux at the black brane horizon, and the nonextremality 
parameter (temperature). Moreover, for the correct thermodynamics, the five-form 
flux at the horizon must depend on the temperature. This was explicitly demonstrated 
in [21]. The correct horizon flux/temperature relation for the cascading gauge theories 
was determined without holographic renormalization in [20] ^'^ by requiring that the 
glueball mass scale of the theory is kept fixed. The same conclusion could be reached^^ 
by keeping the 'holographic scale' (the radial coordinate) 'fixed' — independent of the 
nonextremality parameter. One way to achieve this is to require that the asymptotic 
relation between two supergravity modes (which are affected by the temperature only 
at the subleading order near the boundary) are kept fixed — the same as in the extremal 
(super symmetric) case. In the case of nonextremal MN deformation, this would be the 
statement that the asymptotic relation between the 5*2 radius and the dilaton near the 
boundary of (jB)) is the same as in the extremal MN solution. For large oq, this leads to 
fl33|) . In our analysis of the MN thermodynamics we assume ()33|) to be correct. Then, 
using dnni), (EH), (El we find 

2 _ - ^ ^ - ^^^^ - d\nT/dao 
~^ ~ dS ~f dS ~ dlnS ~ dlnS/dao 

= - - + Oia,') . 
"o 

Notice that for ^ 1, is purely imaginary; moreover it vanishes in the limit 
ao —>■ +00, in agreement with the computation for the speed of sound in six dimensional 
LST [30]. In the next section we reproduce ()40|) directly from the computation of the 
lowest quasinormal mode in finite temperature MN geometry. 

We conclude this section by commenting on the previous proposals for the equations 
of state for the finite temperature MN solution. 

■ In [14] we discussed the thermodynamics of MN solution assuming that 

^'^^ 0. (41) 



da, 







""^^The authors then computed the remaining thermodynamic quantities from (|39|l . 
-•^^I would hke to thank Ofer Aharony for discussing this point. 
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We argued here that PTj) is incorrect, rather one should use The speed of sound 
deduced under the assumption (|¥T|l disagrees with the exphcit computation of the 
lowest quasinormal mode dispersion relation of the next section. 

■ In [15] the authors computed the free energy of the finite temperature MN geometry 
using a reference background subtraction prescription [31]. It is known that background 
subtraction as a method for computing the free energy does not work for charged black 
holes in AdS^, and for the supergravity dual to mass deformed J\f = A SYM theory 
[32-34]. We claim here that is does not work for the nonextremal MN background as 
well. The conclusion of the background subtraction method to finite temperature MN 
geometry was (see Figs. 20,21 of [15]) that in the limit S +oo (the limit of standard 
flat NS5 brane thermodynamics) both S oc S and —P = JF oc —S. Thus in the high 
energy (entropy) limit S cc P. So in the limit S —>■ +oo 

dV 

v1 = constant ^ , 

which contradicts the fiat LST result [30]. 

4 Sound wave in nonextremal MN geometry 

The general prescription for the evaluation of the quasinormal spectra of the nonex- 
tremal gravitational backgrounds was explained in [8]. The application of the pro- 
cedure to a particular solution is very straightforward. The discussion here closely 
resembles [22], so we outline only the main steps. 

Consider fluctuations in the background geometry |T7j) 

dfiu ^ 9fj.u ~^ , 

g^g' + gi, (42) 

where {g^^y, g^, $''} are the black brane background configuration (satisfying (pi?|) . (fT!Hl ). 
and {hf^i^,gi, $1} are the fluctuations. We choose the gauge 

htr = hxir = hrr = 0. (43) 

Additionally, we assume that all the fluctuations depend only on (t, 0:3, r), i.e., we have 
an 0(2) rotational symmetry in the X1X2 plane. At a linearized level we find that the 
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following sets of fluctuations to decouple from each other 

{htxi , hxixs} ; 
{htx2J hx2X3} ; 

{htt: haa = hx^xi + h 



(44) 



The last set of fluctuations is a holographic dual to the sound waves in MN model 
which is of interest here. Introduce 



htt 






2 

3e 


4 $6 






htz 


^—iujt+iqx3 


n''\9') 


2 

3e 


3^ 




haa 


^—iidt+iqxs 


n''\9') 


2 

3e 


3 


Haa 1 


hzz 


^-ioJt+iqxs 


n''\9') 


2 

3e 


3* 


Hzz , 


91 


^^-iut+iqx; ^ ^ 












$1 


^^-iut+iqx, ^ ^ 













(45) 



where {Hu, Ht^, Haa, H^z, 4',4'} ci^e functions of a radial coordinate only. Further we 
introduce diffeomorphism invariant fluctuations 



=4— Htz + 2 Hzz — Haa 1 n fi 



UJ 



g2 2fb9b'^l- fb9'b-39bfi\ , 



+ '^^fb Htt 



39b9'b 



^9'b-^9b^'b 
39b^'b 



Ha 
Haa ■ 



(46) 



The quasinormal mode spectrum is determined [8] by imposing the incoming boundary 
condition at the horizon, and the Dirichlet condition at the boundary on the diffeomor- 
phism invariant fluctuations ()46p . So at the horizon (/^ 0+) we have (see also [35,22]) 



Znir) = h'-^ ZH{r) , Z,{r) = f, z,(r) , Z^(r) = /, ^$(r) , (47) 



'2nT 
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where {zh, Zg, z^} are regular at the horizon; while at the boundary (/;, — 1_) 

zh ^ 0, Zg ^0, 2;$ . (4J 
Taking the hydro dynamic limit, 



00^0, g^O, -^Vs, (49) 

Q 

we find the following set of equations (in the radial coordinate ((221)) 
=z'Ij + Ah z'jj + Bh Zh + Ch Zg + Dh z^ , 

=Z'^ + AgZ'g+ Bg Zg ^Cg Z^+ Dg {X Z'^ " Z h) , (SO) 

=z'^ + Aq, z'i^ + Bq, + Zg + Dq, {x z'jj - Zh) , 

where the coefficients are given explicitly in the Appendix. The symmetry of these 
coefficients imply that the equation for 

X = Zg — Zii, (51) 

is particularly simple^^: 

= xX" + X' . (52) 

A unique solution of consistent with boundary conditions (jlZj), (jlHj) is A(x) = 0. 
Thus we conclude that in the hydrodynamic approximation 

Zg = z<j) . (53) 

We now proceed with solving ()50|) . subject to the boundary conditions ()47p . 
in the approximation oq ^ 1. Using results of section 3.2 and (|53|). we find 

+ f l-3^^-2t;g _ 16x{l-vl) 1 
U-^^^ + \3;(_2t;2 + i + a;2) a^(-2t;2 + l + ^T) ^ ° ' i 

32(1 -2f 2) 32(3 + 2(1 - x^) ln(l -x^) -6f2 + x2(l + 2t;2)) 



2 



a2(-2t;2 + 1 + x2)t.2 a^(-l + x2)(-2t;2 + 1 + x^)v 



n // 1 / f 8 16(ln(l - a;2) - 3) 

° + X + 1 ag(l - x2)2 + (1 - xr< ^ 

4:vI{xz'h - zh){x^ + ln(l - x2)(l - x^)) , _2^ 



-6^ 



a^x4(-2t;2 + l + x2)(l-x2) + 

Notice that eqs. are perturbative in Oq^, but exact in Vg 

^^This is true only in the hydrodynamic approximation (|49f) . 



(54) 



(55) 
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4.1 Leading order solution of eqs. (jMj) . (fS^ 

To leading order in a^^ we have from (j55|l 



= 4' + -4, (56) 



which subject to the boundary conditions (jTTj) . (jiH|) has a unique solution 



pleading 
^9 



(57) 



To leading order in ^ we have from ()54p 



= ^/f + ^ 7, -tt: z'jr H zh . (5^ 



The general solution of (j58p takes the form 

ZH = Ci \ l-x^- 2vH + C2 I (1 - - 2t;2) Inx + 2 - 4^^ [. , (59) 



where Ci are the two integration constants. The incoming wave boundary condition 
(|47p implies that C2 = 0. The Dirichlet condition at the boundary then determines 

[vl^^d'^'s^^ = , (60) 

which is consistent with results of [30]. Without loss of generality we can set Ci = 1, 
thus 

4'""'''' = l-x\ (61) 
4.2 Next-to-leading order solution of eqs. (jMj) . (IS^ 

The next-to-leading corrections to ^^adrn3| yavsi satisfy the Dirichlet condition 

both at the horizon and the boundary. 

Analyzing ()55|) it is easy to see that next-to-leading correction to Zg^"''^^^^, z^^^ oc 
^ %' . It satisfies the following equation 



n untiy, ^ 1 .ntiy ^ 4(t;f)^(x^ + ln(l-x^)(l-x^)) 
°=t"^ ] +^ ] + a^x^(l-x^) ' 



where we substituted the leading solution for zh, (jHH), and used the fact that the speed 
of sound vanishes to leading order in Oq^, (lUUj). A unique solution to (jU^ satisfying all 
boundary conditions takes the form 

~ at x^ ■ 
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Next, consider next-to-leading correction to 2;^'^'^*"^. Again, given the leading order 
solutions and (jU^ . it is clear from that z^f- oc Oq^, (t>"*')^ oc Oq^. Setting 



we find from (j54j) 



4 ' 



(64) 



0H.^T + ^[.^^1' + ^.]^*' + ^. (65) 
The general solution to ()65p takes form 

2^*' = |l -a;^| +C4 |(1 - x^)lnx + 2| - 4-2/3, (66) 

where Ci are the two integration constants. Regularity of the solution at the horizon 
implies that C4 = 0; the Dirichlet condition at the boundary implies (3 = —2; finally, 
the horizon boundary condition on the next-to-leading correction implies C3 = 0. Thus 
we conclude^^ 

< = 0{a^') , (67) 

and 



= -^ + 0{a^') , (6^ 



in precise agreement with ()40p . 



5 Conclusion 

In this paper we presented a simple proof of the Gubser-Mitra conjecture [3] for the 
supergravity backgrounds with translationary invariant horizons what can be inter- 
preted as holographic duals to strongly coupled finite temperature gauge theories. We 
illustrated the discussion with the study of thermodynamics and the hydrodynamics 
of the MN [6] geometries. We explicitly demonstrated that finite temperature MN 
background in the near Hagedorn regime has a negative specific heat, and the tachy- 
onic lowest quasinormal mode. Moreover, the specific heat and the speed of sound 
(extracted from the dispersion relation for the lowest quasinormal mode) satisfy (0), 
which is predicted by our proof of the Gubser-Mitra conjecture. 

^^In order to evaluate zj*' explicitly, one needs to know n = 3 functional coefficients of (|27|l . This 
is not necessary for the evaluation of the speed of sound to first nonvanishing order. 
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Since small fluctuations of a thermodynamic system are described by hydrodynam- 
ics, it is an interesting question as to whether there are other instabilities (besides the 
sound mode) in the thermodynamic system with a negative specific heat^^. We do 
not know the answer in general, however, in the specific example of the MN model, 
we believe that there are none. First, the shear viscosity of the MN system is uni- 
versally [36-39] S'/(47r); thus, it is positive. Second, both the bulk viscosity and the 
R-charge diffusion constant of the LST [30] is finite, and positive. Our analysis sug- 
gests that the ag —>■ +00 limit of the MN model smoothly reproduces six-dimensional 
LST results, thus we find it unlikely that the bulk viscosity or the R-charge diffusion 
constant would be negative for large Qq. 
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Appendix 

Coefficients of fluctuation equations (j^Uj) : 

{x{2x%g,-g,-xg',)-vl{-g', + 2%g,))x{-g', + 2%g,) ' ^ > 



^ {-g'^ + 2^',g,){x{2x%g,-g,-xg',)-vl{-g', + 2%g,))' ^ ^ 



would like to thank Andrei Starinets for discussing this. 
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c = I 

" 3gb{-g', + 2^',g,){2x^^',g, - xg, - x^g', + vl{g', - 2%g,)){-l + gl + 2g,)vl 

X {^x{-3xg'^gl - Sxg', + Ux^'^g^ - 12x<^',gl - 8g, + 2gl + lQgl){-m'^gl 
+ Ag,g', + M'^gl - Sxg',^',g, + xg',^) + vl{-Z2g'^gl + 512*^52 ^ ge^^^^^f 

- ^Qxg'^^',gl + mx^'^g',gl + 48x5^ ^'b^^' - ^^^^9bK9b - 384x$f ^^^^ 

+ 128^,^$;^^ + ?,Ugtx^'i - bl2^'M + 640$;^;^,^ - ?,Mx^'^gl + I2xg',^gl 

- l2Sgl^'^ - U^gt^'^ + I2xg'^ + 1285;'5(. - 160^^ ^,2)^ , 

(71) 



r, _ 8x(8$;ff,^ - Ag,g', - Sx^f^g + %xg',%g, - xg'^) + 32(-ff; + 2<^[g,fvl 

Z{-g', + 2%g,){2x^%g,-xg,-x^g', + vM-2%g,))vl ' ^ ^ 

^ = = - , (73) 



^ -8$;^," + Ag^g[ + 8x$f - 8xg;$;gfe + xg'^ 

9 * Qxgl{-1 + ^2 + 2^,)(2x2$;^, - xg, - x^t, + v'M - '^%9b)){-9i + '^%9b) 

X (x{24x%^gl - 48a;$f + ^'^xg'.'^'.g, + 32<^',gl - 20<l>',gl + 3xg',^gl - 3xg'^ 



- 8g,g'b - H9I) + i^^^Ul - 24$f - 24$^^^^, + - H^gDv] , 



C„ —Ba 



(74) 

{-S^'^gl + 4(7b(7^ + M'jgl - Sxg'^^'^g^ + xg';}){A^[gl - A%gi, + + g'^) 
3g,{-gi + 2<^',g,){2x^<^>',g, - xg, - x^, + vl{g', - 2%g,)){-l + gl + 2g,) ' 

(75) 



8x2^;.(-l + ^2 + 2g,){2x^%g, - xg, - xY, + - '^%9b)){-9b + '^%9b) ' 

(76) 



19 



References 



[1] J. M. Maldacena, "The large N limit of superconformal field theories and super- 
gravity," Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)] |arXiv:hep-th797Tl200j . 

[2] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large 
N field theories, string theory and gravity," Phys. Rept. 323, 183 (2000) 
|arXiv:hep-t h/99051ir|. 

[3] S. S. Gubser and I. Mitra, "Instability of charged black holes in anti-de Sitter 
space," .arXiv:hep-th/0009126| 

[4] H. S. Reall, "Classical and thermodynamic stability of black branes," Phys. Rev. 
D 64, 044005 (2001) [arXiv:hep-th/0104071 . 



[5] S. S. Gubser and A. Ozakin, "Universality classes for horizon instabilities," JHEP 
0305, 010 (2003) arXiv:hep-th/0301002^ . 

[6] J. M. Maldacena and C. Nunez, "Towards the large N limit of pure N = 1 super 
Yang Mills," Phys. Rev. Lett. 86, 588 (2001) |arXiv:hep-th/000800l|. 

[7] O. Aharony, S. Minwalla and T. Wiseman, "Plasma-balls in large N gauge theories 
and localized black holes," arXiv: hep-th / 05072 19| 

[8] P. K. Kovtun and A. O. Starinets, "Quasinormal modes and holography," 
|arXiv:hep-th/0506184| 

[9] M. Berkooz, M. Rozali and N. Seiberg, "Matrix description of M theory on T**4 
and T**5," Phys. Lett. B 408, 105 (1997) |arXiv:hep-th/ 9704089|. 

[10] N. Seiberg, "New theories in six dimensions and matrix description of M-theory 
on T**5 and T**5/Z(2)," Phys. Lett. B 408, 98 (1997) | arXiv:hep-th/9705221J . 

[11] A. H. Chamseddine and M. S. Volkov, "Non-Abelian BPS monopoles in N = 4 
gauged supergravity," Phys. Rev. Lett. 79, 3343 (1997) |arXiv:hep-th/9707176| . 

[12] A. H. Chamseddine and M. S. Volkov, "Non-Abelian solitons in N = 4 gauged 
supergravity and leading order string theory," Phys. Rev. D 57, 6242 (1998) 
|arXiv:hep-th/9711181 . 



20 



[13] A. Buchel and A. R. Frey, "Comments on supergravity dual of pure N = 1 super 
Yang Mills theory with unbroken chiral symmetry," Phys. Rev. D 64, 064007 
(2001) Miv:hep-th/0103022| . 

[14] A. Buchel, "On the thermodynamic instability of LST," ^arXiv:hep-th/0107102 

[15] S. S. Gubser, A. A. Tseythn and M. S. Volkov, "Non-Abelian 4-d black 
holes, wrapped 5-branes, and their dual descriptions," JHEP 0109, 017 (2001) 
)arXiv:hep-th/0 108205 1 . 

[16] I. R. Klebanov and A. A. Tseytlin, "Gravity duals of supersymmetric SU(N) x 
SU(N+M) gauge theories," Nucl. Phys. B 578, 123 (2000) arXiv:hep-th/0002159j . 

[17] I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: 
Duality cascades and chiSB-resolution of naked singularities," JHEP 0008, 052 
(2000) |arxiV:hep-th/0007T9ll . 

[18] A. Buchel, "Finite temperature resolution of the Klebanov-Tseytlin singularity," 
Nucl. Phys. B 600, 219 (2001) |arXiv:hep-th/0011146 . 

[19] A. Buchel, C. P. Herzog, I. R. Klebanov, L. A. Pando Zayas and A. A. Tseytlin, 
"Non-extremal gravity duals for fractional D3-branes on the conifold," JHEP 
0104, 033 (2001) |arXi v:hep-th /010210j . 

[20] S. S. Gubser, C. P. Herzog, I. R. Klebanov and A. A. Tseytlin, "Restora- 
tion of chiral symmetry: A supergravity perspective," JHEP 0105, 028 (2001) 
|arXiv:hep^h/0102172^ . 

[21] O. Aharony, A. Buchel and A. Yarom, "Holographic renormalization of cascading 
gauge theories," Phys. Rev. D 72, 066003 (2005) |arXiv:hep-th70506002| . 

[22] P. Benincasa, A. Buchel and A. O. Starinets, "Sound waves in strongly coupled 
non-conformal gauge theory plasma," arXiv:hep-th/0507026 

[23] J. M. Maldacena and A. Strominger, "Semiclassical decay of near-extremal five- 
branes," JHEP 9712, 008 (1997) |arXiv:hep-th/9710014 . 

[24] M. Henningson and K. Skenderis, "The holographic Weyl anomaly," JHEP 9807, 
023 (1998) |arXiv:hep-th/98060871. 



21 



[25] V. Balasubramanian and P. Kraus, "A stress tensor for anti-de Sitter gravity," 
Commun. Math. Phys. 208, 413 (1999) |arXiv:hep-th/9902121j. 

[26] S. de Haro, S. N. Solodukhin and K. Skenderis, "Holographic reconstruction 
of spacetime and renormahzation in the AdS/CFT correspondence," Commun. 
Math. Phys. 217, 595 (2001) |arXiv:hep-th/0002230. . 

[27] M. Bianchi, D. Z. Freedman and K. Skenderis, "How to go with an RG flow," 
JHEP 0108, 041 (2001) |arXiv:hep-th/0105276| . 

[28] M. Bianchi, D. Z. Freedman and K. Skenderis, "Holographic renormahzation," 
Nucl. Phys. B 631, 159 (2002) |arXiv:hep-th/0112119| . 

[29] I. Papadimitriou and K. Skenderis, "Thermodynamics of asymptotically locally 
AdS spacetimes," JHEP 0508, 004 (2005) |arXiv:hep-th/0505190| . 

[30] A. Parnachev and A. Starinets, "The silence of the little strings," 
|arXiv:hep-th/0506144| 

[31] S. W. Hawking and G. T. Horowitz, "The Gravitational Hamiltonian, ac- 
tion, entropy and surface terms," Class. Quant. Grav. 13, 1487 (1996) 
|arXiv :gFqc79501014| . 

[32] A. Buchel and L. A. Pando Zayas, "Hagedorn vs. Hawking-Page transition in 
string theory," Phys. Rev. D 68, 066012 (2003) arXi v:hep-th/0305179| . 

[33] A. Buchel and J. T. Liu, "Thermodynamics of the N = 2* flow," JHEP 0311, 031 
(2003) |arXiv:hep-th /0305064 . 

[34] A. Buchel, "N = 2* hydrodynamics," Nucl. Phys. B 708, 451 (2005) 
|arXiv:hep-th/0406200| . 

[35] G. Policastro, D. T. Son and A. O. Starinets, "From AdS/CFT corre- 
spondence to hydrodynamics. II: Sound waves," JHEP 0212, 054 (2002) 
[arXiv:hep-th/0210220j . 

[36] P. Kovtun, D. T. Son and A. O. Starinets, "Holography and hydrodynamics: 
Diffusion on stretched horizons," JHEP 0310, 064 (2003) |arXiv:hep- th/0309213|. 



22 



[37] A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," 
Phys. Rev. Lett. 93, 090602 (2004) |arXiv:hep-th/0311175| . 

[38] P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting 
quantum field theories from black hole physics," Phys. Rev. Lett. 94, 111601 
(2005) |arXiv:hep-th/040523ll . 

[39] A. Buchel, "On universality of stress-energy tensor correlation functions in super- 
gravity," Phys. Lett. B 609, 392 (2005) |arXiv:hep-th /0408095 . 



23 



